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Abstract 
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I. INTRODUCTION 


The standard three-flavor paradigm of neutrino mixing has been solidly established by the solar, 


atmospheric, reactor and accelerator neutrino e^eriments [Ij]. However, the anomalies from the 


LSND 


fl, 


MiniBooNE Q], reactor antineutrino and Gallmm radioactive source [5| experiments 


Q. 


B 


as well as the analysis of current cosmological observations [6| have provided interesting hints that 
the three-flavor framework may be incomplete and there may exist additional light sterile neutrinos 
[ 7 , 8|. From the theoretical point of view, sterile neutrinos are regarded as natural ingredients of 
different types of seesaw models 

number of sterile neutrino species, the neutrino mass scale as well as the patten of active-sterile 
mixing are all currently unknown, which consistent one of the most important tasks for future 
precision neutrino experiments. 

The mixing of three active neutrinos and with the possible existing sterile ones could 


lepton mixing matrix 


ix Q 


result in unitarity violation (UV) Ill4l5l| of the 3x3 Maki-Nakagawa-Sakata-Pontecorvo (MNSP) 


where the sterile neutrinos can be either heavy or light. We call the 


former case indirect unitarity violation (lUV) and the latter direct unitarity violation (DUV) {ll|. 
These two kinds of unitarity violation may have distinct effects in neutrino oscillation experiments, 
where the reason is that heavy^ massive sterile neutrinos are kinematically forbidden in the neutrino 
production and detection processes of low energy neutrino experiments, while light sterile neutrinos 
are able to participate in neutrino oscillations as their active partners. The purpose of this paper is 
to derive a complete set of approximate analytical formulas for neutrino oscillation probabilities in 
matter in the presence of both the light and the heavy sterile neutrinos. The different effects of the 
DUV and the lUV on neutrino oscillation experiments are particularly emphasized. A complete 
numerical calculation of these probabilities is also presented to make a comparative study. 

This paper is organized as follows. In Sec. II, we briefly introduce the general framework of 
the (3-|-l-|-2) scenario, which is essential for studying the UV effects. The different effects of the 
DUV and the lUV are discussed. In Sec. Ill, we present the series expansion formulas for the 
neutrino oscillation probabilities in matter of constant density up to the first order in both the 
mass hierarchy parameter a = Am 2 i/Am'^i and those UV parameters sfj (for i = 1,2,3, and 
j = 4,5,6). The corresponding probabilities in vacuum up to the same order are also presented. 


^ Our definition of “heavy” sterile neutrinos assumes that they cannot be produced via the weak interaction process 
as a result of energy conservation. Therefore, the minimal masses of the heavy neutrinos are process dependent. 
For illustration, in the on-shell decay, W'^ —^ U -|- the maximal mass of neutrinos which can be produced 
is 'Jrnd/ff — rnf. In this particular process the definition of ’’heavy” is > \/'m^ — rnf. 
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Section IV is attributed to the numerical calculations. We show in this section the qualitative 
behaviors of neutrino oscillation probabilities, and test the accuracy of these analytical formulas. 
Finally, we summarise in Sec. V, and the details of the perturbative expansion of the neutrino 
oscillation probabilities in the case of UV are given in the appendices. 


II. GENERAL FRAMEWORK IN THE PRESENCE OF UNITARITY VIOLATION 

In this paper, we consider the (3+1+2) scenario with three additional species of sterile neu¬ 
trinos, of which one is the light sterile neutrino and the other two are heavy sterile ones and 
discuss the possibility of determining both the DUV and the lUV parameters in neutrino oscillation 
experiments. 

A. Parametrization of the Lepton Mixing Matrix 

In the (3+1+2) scenario, the full picture of neutrino mixing are described by a 6 x 6 unitary 
matrix which can be decomposed as [iT, ^ 


U = 


1 0 

0 Vr 



AVn R 


.U,SV, U,B. 


( 1 ) 


where 1 + and Uq are unitary matrices while A, B, R and S are not unitary, 0 and 1 stand for the 


zero and identity matrices respective 
these matrices can be found in Ref. 


y. All of them are 3x3 matricies. The explicit expressions of 


171 ] . The full 6 x 6 unitary mixing matrix U is parameterized 


by altogether 15 mixing angles and 15 CP-violating phases (including the Majorana phases). 

In the presence of three additional sterile neutrinos, the Lagrangian of the standard weak 
charged-current interaction should be written as 


-r =A. 

V2 


e ^ 1L 







V 

^2 

+ R 

^5 




L 

W 

L- 


W- + h.c. , 


( 2 ) 


where 1^4 is the mass eigenstate of the light sterile neutrino and 1 ^ 5 , are mass eigenstates of the 
two heavy sterile neutrinos. V = AVq is just the MNSP matrix, which is in general non-unitary 
in this scenario. Equation (2) indicates that the sterile neutrinos can participate in the charged- 
current interaction through their mixing with the active ones (which is described by the mixing 


Current cosmological observations favor at most one species of light sterile neutrinos 
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matrix R) and therefore result in low energy signals of the UV. The matrix Vq in Eq. (1) can be 
parametrised using the standard parametrization as [l| 

/ 


En = 


Cl oC 


12^3 


^12*^13 


«i3e 


-i<5\ 


o*<5 

^i 8 


\ ^ 12^23 ^12^23^13^^° ^12^23 ^12^23^13^^^ ^23^13 J 


( 3 ) 


To the order of sfj {i = 1, 2,3 and j = 4, 5,6), A and R can be approximately written as 


Q 


T ~ 1 — 


Tt (®14 + '5 i 5 + 'Sie) 


Sl4'S24 + SiR'S'jR + SIR'S 


15'’25 


’16'’26 


(s24 + S 25 + 'sie) 


0 

0 


\ 


Si4'S34 + S 45 S 35 + S^gSgg S24'S34 + S25'S35 + S26'S36 


-*- / 2 , 2 , 2 \ 

- (S34 + S35 + S36) y 


( 4 ) 



fs* 

.Si4 

s* 

*16 

R ~ 

«24 

*25 

*26 



*35 

*36y 

Here = sinO^j 


= C 


( 5 ) 












phase angles, respectively. We can find that there are altogether 9 additional mixing angles and 
9 additional phases in R and A, which are relevant to the low energy experiments. These UV 
parameters can be divided into two groups: 

• (^ 14 ) ^ 24 ) ^ 34 ) ^^*4 (<^ 14 ) <^24 *^ 34 ) parameters of the DUV describing the mixing between 
three active neutrinos and the light sterile neutrino. Since the light sterile neutrino can be 
produced and detected in the low energy experiments, these parameters will also appear in 
the neutrino oscillation probabilities as those standard mixing parameters. 

• (^ 15 > ^ 16 > ^ 25 > ^ 26 > ^ 35 > ^nd <^ 25 > <^ 26 > ^ 35 > <536) are the lUV parameters that 

describe the mixing between three active neutrinos and the heavy sterile ones. The lUV 
affects the processes of neutrino production and detection, and therefore induce the “zero- 
distance” effects [l^ in neutrino oscillations. 

The analyses in Refs, [l^, [l^ show that leptonic and semileptonic decays are quite sensitive to the 
lUV of the leptonic mixing matrix. The strongest constraints come from the rare charged lepton 
decay experiments, which indicate that the combinations < 1%, with a,l3 = e,/r, r, 

should be satisfied On the other hand, the DUV parameters can be revealed in the very-short- 


® Note that, instead of the 3x3 non-unitary MNSP matrix E, the matrix U is the 3x4 left-up sub-matrix of W as 
showed in Eq. (6). Since the light sterile neutrino can be kinematically produced in the electroweak decays, the 


DUV parameters are irrelevant to the strong constraints obtained there [h IQ. 
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baseline neutrino oscillation experiments. The global analysis of current results of the short-baseline 
neutrino experiments points towards the eV scale feht sterile neutrino(s) with active-sterile mixing 
parameters |C/q, 4 P (a = e,/i,r) of a few percent [Sj]. In this paper, we simply take all the active- 
sterile mixing angles {i = 1,2,3 and j = 4,5, 6 ) as small parameters of the same order s. 


B. Neutrino Oscillation in the Presence of Unitarity Violation 


In this paper, we restrict us to the typical neutrino oscillation process ^ where both 
the production of and the detection of are via the charged-current interaction. Since in 
our proposed (3-|-1-1-2) scenario, there is only one light sterile neutrino that can participate in 
the neutrino oscillation together with the three active ones, only the elements in the 3x4 left-up 
sub-matrix U of the full 6 x 6 mixing matrix Li are related to the neutrino oscillation probabilities: 




VrJ 


= u 


^2 

^3 


/ c /,1 c /,2 

U^l U^2 U^3 U^4 
C/,2 U,3 






In this case, the neutrino oscillation probability in vacuum should be written as [l^, ^ 


P ( ^ ^ p) — 


(^U*e-iELuT) 


q/3 


E 

i=l,2,3,4 




0.(3 


(C/C/t)„„ (c/c/t) 


/3/3 


E 

*=1,2,3,4 


E E ic'/sii 

*=1,2,3,4 *=1,2,3,4 


E 

*=1,2,3,4 

2 


*=1,2,3,4 


-4 E Re Sin 2 A^, ± 2 E Im sin 2 A^, 

j>i j>i 

2 

1 


E ClVa, 

E ( E lEii^) ^ 

*=1,2,3,4 *=1,2,3,4 

-4Re sin2 A 21 ± 2Im [U sin 2 A 2 i 

-4Re [C/„iC/^ 3 C/* 3 C/^i] sin2 A 33 ± 2Im [U sin 2 A 34 
-4Re [C/„2C/;33t^:3t^;32] sm' A 32 ± 2Im [U sin2A32 


( 6 ) 
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-4Re sin2 A 41 ± 2Im sm2A,, 

-4Re [U^2Uf,,U*4Uk] sin^ A 42 ± 2Im [U^2U^,U*4U*p2] 

-4Re [U^sU^,K,U*^,] sin2 A 43 ± 2Im sin 2 A 43 | , (7) 

where Aj^ = Am'j^L/AE with Am^j = — mf being the neutrino mass-squared difference. Here 

the Greek letters a, (3 are the flavor indices e, //, r, while the Latin letters i, j are the indices 
of mass eigenstates. Note that the indices i, j run over only the light neutrinos (both the active 
and sterile ones) which can be kinematically produced in neutrino oscillation experiments, and 
the normalization factor l/(X]j=i 2 3 4 iKiil^)(Z]i=i 2 3 4 ensures that at the source we have 

| 44 -(VL''' —)• = 1 and at the detector \A{i 2 -W~ = 1 is satisfied. 

When neutrinos propagate through the Earth before reaching the detector, neutrinos may inter¬ 
act with electrons, protons and neutrons in the medium via the charged-current (CC) interactions 
or neutral-current (NC) interactions. The coherent forward scattering from the constituents of the 
Earth matter will modify the evolution behavior of neutrinos. For the propagation of neutrinos in 
matter, the Hamiltonian can be written in the basis of mass eigenstates in vacuum as 

n = E + U^AU* , (8) 

where E = diag {E 4 , E 2 , £' 3 , £ 4 } is the energy matrix, A = diag {V^c ~ ~^nc^ ~^nc} 

= 'Y^GpTi^ and = Gpn^lV2 (n^ and are the electron and neutron number densities, 
respectively) are the CC and the NC contributions to the neutrino matter potentials respectvely. 
U is just the 3x4 non-unitary mixing matrix given in Eq. ( 6 ). 

Note that in Eq. ( 8 ), £ is a 4 x 4 diagonal matrix because altogether four mass eigenstates i^i, 
1 ^ 2 , z ^3 and 2/4 can be kinematically produced and therefore participate in the neutrino oscillation, 
while A is a 3 X 3 matrix since only three active left-handed neutrinos zZg, and are involved 
in the CC or the NC interactions. The Hamiltonian in Eq. ( 8 ) holds for neutrinos, whereas one 
has to perform the replacements U ^ U*, V(jp, ~^cc '^NC ~^NC antineutrinos. 

In the case of matter of constant density, the Hermitian matrix T-L can be diagonalized by 
a unitary transformation H = XEX\ where £ = diag |£ 4 , £ 2 , £ 3 , £ 4 ! is the effective energy 
matrix in matter and X is a 4 x 4 unitary matrix. It can be inferred that the matrix U = UX*, 
which is also a 3 x 4 non-unitary matrix, can be regarded as the effective leptonic mixing matrix 
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in matter. Therefore, we can write down the neutrino oscillation probabilities in matter as 

2 


P ( ^ a~^ ^ p) — 


1 


E E 


i=l,2,3,4 


i=l,2,3,4 


i=l,2,3,4 


-4ERe 

j>i 


UaiUf^j%U*pi 


sin^ Ajj ± 2 E I™ 

j>i 


Uo^iUpjKMi 


sin2A i \ ,(9) 


where /AE with Am^^ = — m? = 2E^{E- — E^) is the effective neutrino mass- 

squared difference in matter. 

One can find from Eqs. (7) and (9), the UV in the MNSP matrix may result in three kinds of 
different effects on the neutrino oscillation probabilities summarised as follows. 


The lUV can affect the neutrino production and detection, and then generate the “zero- 
distance” effect 1^, i.e., at L = 0 it follows that 

2 


P p)\l=0 = 




{uux. 


^ 0 (a / /3) , 


( 10 ) 


/3/3 


where U is the 3x4 matrix in Eq. (6). It means the flavor transition already took place at 
the source before the oscillation begins, with very tiny transition probability of the order 
It is worth to mention that the “zero-distance” effect is irrelevant to the DUV, which means 
this effect would not be observed if there exists only the light sterile neutrinos. 


• The existence of light sterile neutrinos will introduce both additional CP-conserving and CP- 
violating oscillatory terms in neutrino oscillation probabilities. Their oscillation amplitudes 
are of the order and the frequencies are proportional to the newly introduced mass-squared 
differences. The oscillatory behaviors of these new frequencies may be observed at very-short 
baseline neutrino oscillation experiments [^. 


• In the propagation, both the DUV and the lUV can modify the amplitudes of the standard 
oscillatory terms as well as affect the matter effects. Moreover, for the neutrino beams with 
relative hight energies, the non-unitary effects may be largely enhanced by the neutrino- 
matter interactions [2l[ |. 

Although by using Eqs. (7) - (9), the oscillation probabilities can be numerically calculated with 
any required accuracy, the analytical approximate formulas can be useful in showing the UV effects 
in a more transparent way. In the next two sections, we are going to discuss both the analytical 
expansions and numerical calculations of the neutrino oscillation probabilities in the case of UV. 
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Indeed, there are several previous publications discussing the hybrid scenarios with both light 
and heavy sterile neutrinos [^ . however, our study is different because the present work 
presents the long-baseline neutrino oscillation behaviors in the presence of Earth matter effects, 
while the previous ones focus on the property of short baseline neutrino oscillations and the expla¬ 
nation of experimental anomalies. For the first time, we derive a complete set of series expansion 
formulas for neutrino oscillation probabilities in matter of constant density in the presence of both 
light and heavy sterile neutrinos. Formulas for the (3-I-1-I-2) scenario can be easily generalized to 
scenarios with arbitrary numbers of light and heavy sterile neutrinos. Note that sterile neutrino 
induced CP violation may have significant observable effects in current and future long-baseline 


neutrino experiments (e.g., T2K 


M 


251], DUNE 


B), 


which makes our derivations important and 


timely to understand these possible distinct phenomena in long-baseline experiments. 


III. ANALYTICAL EXPANSIONS OF NEUTRINO OSCILLATION PROBABILITIES 


In this section, we derive the series expansion formulas for neutrino oscillation probabilities in 
matter of constant density by using the perturbation theory. The formulas are expanded up to 
the hrst order of both the mass hierarchy parameter a = and the UV parameters 

sfj (i = 1,2,3 and j = 4,5,6). The details of the diagonalization of the Hamiltonian 7{ can be 
found in Appendix A, and the resulting approximate neutrino oscillation probabilities in matter are 
presented in Appendix B, where transition probability from a neutrino 

flavor a to a neutrino flavor /3 in matter. Equations (B3) - (B 8 ) and Eqs. (BIO) - (BI2) can be 
further simplified if one can neglect those terms of 0{si^a), O(sf^sfj) or O(sf^), the approximate 
neutrino oscillation probabilities can then be obtained as: 


P.. 


P, 


MM 


1 - 2 sf 4 - 


1 — 2s24 ~ 


4s 


13 


(1 - ^ccY 


(■Sig + ‘^^b^Nc) sin^(l — ^cc)^3i ; 


( 11 ) 


o „2 1 a(l + Ap^)sin2042sin4023Si3Cos5 
sin Ztf23 U — ^' 524 ; 2 

^CC 


R.A 


+2sin4g23( 


A 


'^a^NC 


CC 


sin^ A 


31 


+ 


‘1®23®13 f ®23®13 Q!Sin2042C23 ,00 1 ^-^a^AfC'^23 , ‘^^b^NC^23 

-h 2 Aio + 


1 “ ^CC V1 “ ^cc 


+ 


4 ' LZ ' 4 14 

^cc ^cc A 

■ [c 23 sin(l -I- A(j(j)/\.^i — S 23 sin(l — sin(l — 

^Cc)^31 


















p. 


efi 


p.. 


p, 


flT 


+A 3 J sin^ 2f?23 f a<a 2 ^ y 


2 I _ 

J~l./^/~1 J. J~\./^/~1 ^ 


^CC 


^cc 


( 12 ) 


1 — 2 S 34 — 


-t“2 sin 4^23 -^23 


+ 


^^ o„2 ^ a(l+ ylp^) sin 2042 sin 4023^13 COS (5 

sin ^(723 — ^^34^ 2 

^a^NC^lS 


A, 


cc 


A, 


+ T^A 


a^NC 


CC 


sin^ A 


31 


4C23S13 f C23S13 , a sin 2042^23 (((^ NC^23 I 2-Rf,2ljY^C23 


1 “ Acc V 1 “ ^cc 


+ 


A. 


A. 


+ 


^cc ^cc 

■ [s^3 sin(l + Acc)^31 - ci3 sm(l - ^cc)^3i] sm(l - ^CC ')^31 


1-^4 


+A34 sin 2023 (^«C 42 - ^ ^ 

4^23^13 / S 23 


A 


CC 


RbANC-^13 _ ^ ) g.^ 2 A 34 , 

1 - A^(j 2 ^ ' 31 ’ 


^ “ Acc I ^ “ ^cc 


(S43 + 2 RfjAj^fj) sin(l — 


+i?42 sin(l — A(^(j)A'^i + 142 cos(l — yl(^(^)A34 
+ ^^3 [(a sin 2042 cos 5 — 2 R^Ap^(j) cos A34 


A 


CC 


+ (a sin2042 sin (5 + 2 Ig.Aj^(j) sin A34] sin(2l(^(^A34)} sin(l — Aq(j)A^i , 


-23 


^^23^13 

1 “ Acc 11 “ Aq(j 


(S43 + 2 Rj^Aj^Q) sin(l — A(jq)A^i 


+i?43 sin(l — A(^q)A^i + 143 cos(l — A34 


’23 


A 


[(a sin 2042 cos S — 2 R^Ajyc) cos A34 


cc 


+ (a sin 2042 sin S + 2 I^A^c) s™ ^31] sin( 2 lcc'^ 3 i)} sin(l - ^c'c)^ 3 i > 

„=„2n/i „2 .2 \ a (1 + ^cc) sin 2042 sin 4023^43 cos 5 

sin Ztl 23 i,-*- “ S24 — S34J 

RnA^jr’S-l 


cc 


+2 sin 4023 ( ^23 + 


''a^NC’^13 I rp A 

- 1 -r ^a^NC 

^CC 


sin^ A 


31 


23 13 54442^2 — A(jq)A^i + 723 sin 2023 sin2A34 


(1 - ^ccY 


+A34 sin^ 20 . 


+ 


23 


2 sin 2023S43 


—ac42 + 


^C'cSl3 I ^bANC^13 


+ 


1 - 1 - 2 I, 


c^NC 


’CC 


cc 


^31 


1 -A 


cc 


a sin 2042 cos 2023 cos 6 


A 


cc 


— 2R^ COs( 2 l(^(^A 34 ) 


cc 


( 13 ) 


( 14 ) 


(15) 
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+ 


a sin 20^2 sin 5 


A 


cc 


— 21^ sin(^^(^A3]^) 


sinA33sin(l - A(jc)^:ii 


~'^^NC sin 26*23533 


fl.COs 2 e ,3 , 


^ “ ^CC V ^cc ^ ~ ^cc 


L 


A, 


■sm(Acc^3i) 


cc 


sinA33sin(l -^^- 0)^31 


( 16 ) 


Here the expressions of i?32) ^13-: ^23^ A2> Asj ^23^ '^a can be found in 

Eqs. (A 2 ) and (A 3 ). Due to the existence of UV (both the direct and the indirect), the “to all” 
transition probability of the a flavor to all active flavors P(y^ —)■ t'e^^j^r) = ^ae + + Par 

a = e,fM,T) is not the unity: 

4si 


P(t'g —)■ 1^6,34,r) ~ 1 — 2S34 + — 


’13 


-A 


[Pb sin(l — A(jq)1^3i 


cc 


+ 4 cos(l — A(j(j)/S. 3 i] sin(l — A^(j)/S. 3 i 


( 17 ) 


P{l 2 ^ b'e^fi^r) ~ 1 — 2^24 + [sin^ 2^23 (s24 “ S 34 ) +-6^23 S™ 4 ^ 23 ] sin^ A 34 +/23 sin 2023 sin 2 A 34 

1 - {-^12 [^*^23 sin(l + A(jfj)IS .33 — §23 sin(l — A(J(j)IS. 3 ^ 

^ ^cc 

~h 2 [^23 COs(l + A(j(j)/S .33 + S 23 COs(l — A(j(j)Ax 3 ^ 

— sin 2023 [-^13 cos(A(^(^A 34 ) + 743 sin(A(^(^A 34 )] sin A 34 } sin(l — A(^(^)A 34 ,( 18 ) 
P{^r ’^e,ix,T) ~ 1 - 2534 - [sin^ 2023 {^24 “ S 34 ) + -^23 sin4023] sin^ A 34 - I 23 sin 2023 sin 2 A 34 

4^23^13 


1 - A, 


CC 


{sin 2023 [-^12 cos(Ac.cA3i) + I12 sin(Ac.c'^3i)] sin A34 


“-^13 [^23 sin(l + A(j(^)/S.3i — C23 sin(l — Auq)/S.3]\ 

+-^13 ['^23 COs(l + A(j(j)/S.33 + C23 COs(l — A(^(^)A34] } sin(l — A(j(j)/S.33 . ( 19 ) 

It is interesting to find that these summed probabilities are all independent of the NC potential 
Vj^(j in these first-order expansions. 

Up to the same order, the corresponding probabilities in vacuum and P{Va ^e,^i,T) i^ 
case of UV are given for comparison: 


-Pee ~ 1 - 2 s ?4 - sin^ 2043 sin^ A 34 , 


( 20 ) 


PfxfM ~ i- 2 s 24 (1 2S24) sin^ 2023 — 4 COS 2023S23S43 — 8P42 COS 2023S23S13] sin^ A34 
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+aA3i sin 2023 ('^12 sin 2023 “ ^ sin 2032S23S13 cos 5 ) sin 2A3^ , 


( 21 ) 


Ptt ~ 1 — 25^4 — [(1 — 2SI4) sin^ 2023 + 4 cos 2023023^^3 

+8 (-Ri2'Si3 + -R 23 S 23 ) COS 2023 C 23 ] sin^ A 31 

+aA3^ sin 2023 (ci2 sin 2023 + ^ sin 2042CI3S43 cos 5) sin 2A3^ , (22) 

Pe^l ~ [s23 sin^ 2043 + 4-^12523^13 + 2aA34 sin 2042 sin 2023S43 sin 6] sin^ A34 

+ (2/42S23S43 + aA34 sin 2042 sin 2023S43 cos 5 ) sin 2A34 , ( 23 ) 

Per ~ [c 23 s™^ 2043 + 4 -Ri 3 C 235 i 3 “ 2aA34 sin 2042 sin 2023^43 sin 6 ] sin^ A34 

+ (2/43C23S43 — 0!A34 sin 2042 sin 2023S43 cos 5 ) sin 2A34 , ( 24 ) 

P^iT ~ sin 2023 { [sin 2023C43 — (524 + SI4) sin 2023 “ 4i?(,S43 
+27223 cos 2023 + 2 q;A 34 sin 2042S43 sin 5 ] sin^ A34 

+ [/23 — aA34 (C42 sin2023 + sin2042 cos 2023S43 cos (5)] sin 2A34I . (25) 

And 


P{^e ^e,a,T) ~ 1 “ 2 si 4 + 4 i 2 feS 43 sm^ A34 + 2^543 sin2A34 , 


( 26 ) 


P{^I1 ~ 1 — 2S24 + [(524 “ 534) sin^ 2023 + 4 (-R12 cos 2023 “ -^13 sin 2023) 523^13 

+7223 sin 4023 ] sin^ A34 + sin 2023 “ 2 -?'i 2 S 235 i 3 ) sin 2A34 , ( 27 ) 

P{^r ~ 1 - 2si4 - [(5^4 - sl^) sin^ 2023 + 4 (-R12 sin 29 ^^ + 7243 cos 2033) C23S43 

+7223 sin 4023 ] sin^ A34 - sin 2023 + 2-?'i3C23Si3) sin 2A34 . ( 28 ) 

The approximate oscillation probabilities in the standard three-flavor case can be easily obtained 
by switching off the active-sterile mixing (i.e., setting all 0 ^^ = 0 , where i = 1 , 2,3 and j = 4 , 5 , 6). 
In this limit, our results are consistent with those obtained in Ref. 271 - 131 ]. 
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FIG. 1: Contours of the probability calculated numerically without any approximations in the standard 
(3+0+0) scenario (top-left) and its corrections by different UV effects — pstandard |.]^g (3_|_0-(_2) 
(bottom-left), (3+1+0) (top-right) and (3+1+2) (bottom-right) scenarios as functions of the neutrino energy 
E and the baseline length L. The probabilities are averaged over a Gaussian energy resolution of 2%. 


IV. NUMERICAL CALCULATIONS AND DISCUSSION 


In this section, we carried out the numerical calculations to discuss: i) the accuracy of above 
approximate formulas and ii) effects of the DUV and the lUV. Here matter of constant density 
3g/cm^ is assumed and the electron fraction is assumed to be 0.5 3^. In our calculation, = 

7.50 X 10-5 eV^ Amia = 2.457 x lO-^ eV^, 0^^ = 33.48°, 0^^ = 8.50°, 023 = 42.3° and S = 306° 
(or Am|i = 7.50 x lO'^ eV^, Am |2 = 2.449 x lO-^ eV^, 0^2 = 33.48°, 0^3 = 8.51°, 023 = 49.5° 
and S = 254°) have been taken for the normal hierarchy (NH) case (or the inverted hierarchy (IH) 


case) 

<^14 = = <5if: = 0°, SoA = = Soa = 60° <5^4 = (5qp; = = 120° huve been chosen as the default 


33l |. For those UV parameters, 0^^ = @24 = ^34 = 7°, 0^^ = @25 = ^35 = ^i6 ~ ^26 ~ ^36 ~ 3°, 

15 = *^16 = *^24 = *^25 = *^26 = ^0° <534 = <535 = dgg 
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FIG. 2: The probabilities at the baseline L = 1000 km calculated numerically without any approxima¬ 
tions as functions of the neutrino energy E in the four different scenarios. The probabilities are averaged 
over a Gaussian energy resolution of 2%. 



FIG. 3: The probabilities {Pep) ^-t the baseline L = 1000 km calculated numerically without any 
approximations as functions of the neutrino energy E in the (3-I-0+0) scenario (solid lines) and the (3+1-I-2) 
scenario (dash-dotted lines) for (a) —>■ oscillation in the NH case, (b) —?> oscillation in the IH 

case, (c) —>• oscillation in the NH case and (d) —>■ oscillation in the IH case. The probabilities 

are averaged over a Gaussian energy resolution of 2%. 
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FIG. 4: The probabilities at the baseline L = 1000 km calculated numerically without any approxima¬ 
tions as functions of the neutrino energy E in the (3-I-0-I-0) scenario (case (a)) and the (3-I-1-I-2) scenario 
with 0^4 = @24 = ^34 = 5° (case (b)), 7° (case (c)) or 10° (case (d)). The probabilities are averaged over a 
Gaussian energy resolution of 2%. 



FIG. 5: The probabilities at the baseline L = 1000 km calculated numerically without any approxima¬ 
tions as functions of the neutrino energy E in the (3-I-0-1-0) scenario (case (a)) and the (3-I-1-I-2) scenario 
with i5i4_i 5 le = 0 , >524 25,26 “ ^34,35,36 ~ 120 ° (case (b)); = 0 , 524 = 525,26 = —60°, 

^34 = 120 ° and 535,36 = “ 120 ° (case (c)); or 5 i4,i5,i6 = 0 , 534,25,26 = “60° and 034 , 35,36 = - 120 ° (case (d)). 
The probabilities are averaged over a Gaussian energy resolution of 2%. 
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FIG. 6: Contours of the probabilities P{y^ fj. t) calculated numerically without any approximations as 
functions of the neutrino energy E and the baseline length L in the (3+0+0) (top-left), (3+0+2) (bottom- 
left), (3+1+0) (top-right) and (3+1+2) (bottom-right) scenarios. The probabilities are averaged over a 
Gaussian energy resolution of 2%. 


values unless otherwise specified. The newly introduced mass-squared difference is assumed 

to be 1.0 eV^, while the masses of two heavy sterile neutrinos are irrelevant to the calculations. 
This particular choice of the UV parameters in our numerical illustration are motivated from 
two aspects. On the one hand, the anomalies in short baseline neutrino oscillation experiments 
snggest that there may exist one or more light sterile nentrinos with the active-sterile mixing of 


0 ( 0 . 1 ) 


7, 


. On the other hand, the mixing of heavy sterile neutrinos with the li ght ones are 


strongly constrained by the electroweak precision data to be less than 0(0.01) Refs. 


13 , 


1^. As 


one can easily image, the larger the active-sterile mixing is, the more significant nnitarity violation 
effects there wonld be. This is why in our (3+1+2) scenario, the DUV effects dominate over the 
lUV effects. 
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FIG. 7: The probabilities P{v^ —>■ ^he baseline L = 1000 km calculated numerically without any 

approximations as functions of the neutrino energy E in the (3+0+0) (solid lines), (3+0+2) (dashed lines), 
(3+1+0) (dotted lines) and (3+1+2) (dash-dotted lines) scenarios for both the NH and the IH cases. The 
probabilities are averaged over a Gaussian energy resolution of 2%. 



E (GeV) 


FIG. 8 : The probabilities P{v^ —>■ at the baseline L = 1000 km calculated numerically without any 

approximations as functions of the neutrino energy E in the (3+1+2) scenario with 9i^ = 024 = ^34 = 5° 
(case (a)), 7° (case (b)) or 10° (case (c)). The probabilities are averaged over a Gaussian energy resolution 
of 2 %. 
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FIG. 9: The probabilities P{v^ ht) s-t the baseline L = 1000 km calculated numerically without any 
approximations as functions of the neutrino energy E in the (3+0+2) scenario with = 0, <525 26 = ^0° 
and (535 36 = 120° (case (1)); or <535 = 0i <525^26 = —60° and <535 gg = —120° (case (2)); the (3+1+0) scenario 

with <534 = 0, 624 = 60° and <534 = 120° (case (3)); or <544 = 0, <524 = —60° and ^34 = —120° (case (4)); 
and the (3+1+2) scenario with <544 45 43 = 0, <524 25 26 — 60° and <534 35 gg = 120° (case (5)); <544 45 4 g = 0, 
<524 = 60°, < 525^26 = —60°, <534 = 120 ° and <535 35 = — 120 ° (case ( 6 )); or <544 45 45 = 0 , ^ 24 , 25,26 = —60° and 
6*34 35 36 = —120° (case (7)). The probabilities are averaged over a Gaussian energy resolution of 2%. 


Taking the transition probability as an example, we illustrate in Fig. 1 the contour structures 
of this transition probability in the standard (3+0+0) scenario (top-left), and its corrections by 
different UV effects (i.e., P^/^ — in the (3+0+2) (bottom-left), (3+1+0) (top-right) and 

(3+1+2) (bottom-right) scenarios as functions of the neutrino energy E and the baseline length 
L. All the probabilities are calculated numerically without any approximations and are averaged 
over a Gaussian energy resolution of 2%. Figure 2 shows how the probability varies with 
the neutrino energy E at a fixed baseline L = 1000 km. We can find that the UV effects could 
modify the oscillation probabilities up to a few percentages while the oscillatory behaviors are not 
largely changed. These features can be understood with the help of Eqs. (14) and (B 6 ), where for 
the long-baseline accelerator neutrino oscillation experiments the fast oscillating terms of A 44 are 
averaged out, and both the DUV and the lUV effects could lead to corrections to the magnitudes of 
the oscillatory terms and small shifts of the oscillatory frequencies (if the matter effects are taken 
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FIG. 10: Contours of the absolute errors of the analytical formulas of in Eqs. (11) - (16) as functions 
of the neutrino energy E and the baseline length L. 

into account) which are all of the order of 

At the baseline L = 1000 km, different oscillatory behaviors of the neutrino and antineutrino 
channels in either the NH or the IH cases are illustrated in Fig. 3. The distinct behaviors between 
the two kinds of mass hierarchies can be observed in the vicinity of the first oscillation maxima (at 
E ^ 2 GeV), and the differences between the neutrino and the antineutrino channels are obvious 
at the second peaks (at E ~ 0.5 GeV). This conclusion is true no matter if the unitarity is violated 
and what kind of UV it is. We study in Figs. 4 and 5 the dependence of the UV corrections to the 
oscillation probabilities on those UV parameters. Generally speaking, the larger the active-sterile 
mixing is, the more significant the UV effects would be. However, whether the oscillatory behaviors 
are slightly enhanced or depressed and the magnitude of the UV corrections depend on the actual 
active-sterile mixing patterns (the mixing angles as well as the phases). 

The equations ^ ^ “ = e, fi, r) always hold if there exists only three 

standard neutrinos and the 3x3 MNSP matrix is unitary. In the presence of UV, these summed 
probabilities may deviate from unity up to the order of s^. The larger the UV is, the more significant 
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the deviation would be (see Fig. 8 ). Moreover, this deviation may be slightly enhanced at the 
region where the matter effects are important. We show in Fig. 6 the contours of P{v^ —)• ^e^ir) 
as functions of the neutrino energy E and baseline length L in the (3+0+0) (top-left), (3+0+2) 
(bottom-left), (3+1+0) (top-right) and (3+1+2) (bottom-right) scenarios respectively. In the DUV 
case (in the presence of light sterile neutrinos), we always have P{v^ < 1. However in 

the lUV case (in the presence of only heavy sterile neutrinos) P{i'^ —)• can be either smaller 

than 1 or slightly larger than 1. This can be easily understood by using Eqs. (17) - (19), where 
the suppression terms (1 — 2 s? 4 ) dominate in the case of DUV. 

To see the intrinsic properties of P{v^ more clearly, we illustrate in Figs. 7-9 the 

one-dimensional plots of this probability as functions of neutrino energy E at the fixed baseline 
of L = 1000 km. Different mass hierarchy scenarios are compared in Fig. 7. One can find the 
curve of the NH case and that of the IH case are approximately horizontally-symmetric, since the 
hierarchy-asymmetric terms dominate the energy dependency as shown in Eqs. (17) - (19). Suppose 
three masses of the active neutrinos are of the normal hierarchy. Fig. 9 shows the probabilities 
P{i^^ —)• ^e,n,T) with different choices of the UV CP-violating phases. In comparison with Fig. 7, 
we can easily find that if the sign of and the signs of all (for i = 1,2, 3 and j = 4, 5,6) 

flipped at the same time, the probabilities P{v^ —)• (fm a = e,/r,T) would approximately 

have the same values. This results can also be easily inferred from Eqs. (26) - (28). 

Before finishing this section, we shall test the accuracy of our approximate analytical formu¬ 
las in Eqs. (11) - (16). We show in Fig. 10 the absolute errors (dehned as the absolute differ¬ 
ences between the analytical and the numerical calculations of the oscillation probabilities, i.e., 
panaiyticai _ pnumericai ^ functions of the neutrino energy E and the baseline length L. We can 
find from the figure that, in the region L/E <C 10^ km/GeV, these analytical formulas are of 
good accuracy (with the absolute errors better than 0.01). Therefore, we conclude that the series 
expansions of the oscillation probabilities where the UV effects have been taken into account are 
useful for understanding of the physics of future long-baseline neutrino oscillation experiments, 
which might be sensitive to the UV effects. 

V. CONCLUSION 

The purpose of this paper is to derive the approximate analytical formulas for neutrino os¬ 
cillation probabilities in matter of constant density in the presence of the direct and the indi¬ 
rect unitarity violation. The probabilities are expanded in both the mass hierarchy parameter 
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a = Am 2 i/Am 3 ]^ and the small active-sterile mixing parameters s'fj (for i = 1, 2,3 and j = 4, 5, 6). 
We have showed that both the charged-current and the neutral-current interactions are non-trivial 
in the presence of UV. Equations (11) - (19) constitute the main analytical conclusions of this 
paper. We have discussed the important properties of the oscillation probabilities by comparing 
among scenarios of different mass hierarchies, different mixing patterns and different CP-violating 
phases. The accuracy of these analytical expressions in different regions of the neutrino energy E 
and the baseline length L are also studied to identify the valid region for realistic applications. Our 
numerical calculations show that, for most of practical long-baseline accelerator experiments, these 
analytical formulas are very good approximations. In particular, we have discussed the deviations 
of the “to all” probabilities P{va from the unity, which provides a definite signal of 

the UV. We find that the behaviors of these summed probabilities are quite different in the case 
of DUV or lUV. However there is an approximate degeneracy if the signs of Aml^ and all the 
UV CP-violating phases flipped at the same time. In the presence of UV, the deviations can be 
enhanced by the matter effects, where only the charged-current interaction is relevant. 

As the measurements of neutrino oscillations enter the precision era, testing the standard three- 
neutrino paradigm becomes a potential probe for new physics beyond the standard neutrino oscil¬ 
lation framework. It is desirable to present a sophisticated study for the test of the UV in future 
long-baseline accelerator neutrino experiments or the Neutrino Factory. The analytical expressions 
given in this paper could reveal the basic dependence of the neutrino oscillation probabilities on 
the DUV or the lUV parameters and help to facilitate the choice of the experimental setup as well 
as the analysis of the data. 

Acknowledgments 

This work is supported by the National Natural Science Foundation of China under Grant Nos. 
11135009, 11305193, and 11105113 , and by the CAS Center for Excellence in Particle Physics 
(CCEPP). 

Appendix A: Diagonalization of the Hamiltonian in matter using the perturbation theory 

In this appendix, we use the perturbation theory to diagonalize the 4x4 effective Hamiltonian 
a. in matter of constant density by the unitary transformation % = XEX"^, where A is a 4 x 4 
unitary matrix and E is diagonal. In the series expansion oiT-L, we regard a = Aml^/Aml^ and 
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s'fj (for i = 1,2,3 and j = 4,5,6) as the small parameters of the same order s^, and perform the 
diagonalization to the first order of a and sjy In the appendices, we shall adopt the following 
abbreviations; 


and 


^ Amqi L 
A 31 = 


AE 


^CC = 




A — 


^11 = 
^22 = 
^33 = 
-^12 = 

-^13 = 
-^23 = 

A 2 = 
As = 
As = 


2 , 2,2 
■Sm + 'S 15 + ^16 ) 

2 , 2,2 

■524 + 'S 25 + ^26 ) 

■534 + 'S35 + -sie ) 

Re (S 14 S 24 + 'Sl5^25 + '^ 16 ^ 26 ) ® 

Re (S 14 S 34 + S 15 S 35 + S 4 gS 3 g) e 

Re [S24S34 + ^25^35 + ^26^36] > 

Im ( 5345*4 + 5455^5 + SigS^g) e-*'^ 
Im (534534 + S 35535 + SigSsg) e 
Im [ 524 S 34 + 525 S 35 + 526 ^ 36 ] ) 


Ai = -^23 cos 26*23 + 2 (^22 “ ^33) 26*23 ) 

A) = -^23 sin 26*23 + '^ 22^23 + As^ls ) 

A = 2 i ?23 sin 2023 + (A2 “ As) cos 2023 ’ 
= “-^12023 + -R13S23 ) 

A = ^2^23 + A 3 C 23 ) 

A = -R12C23 + -^13S23 ) 

A = “A2C23 + A3S23 ) 
h = A 2 S 23 + A 3 C 23 • 


(Al) 


(A2) 


(A3) 


In the basis of neutrino mass eigenstates in vacuum, T-L can be written as 

n = E + U^AU* ~ A • 1 + 77 ^ • diag (O, a, Ami 3 , Aml^) + U^AU* , (A4) 

2 A 

where A = diagjAj^ — RjvC’ ~'^NC^ “Ad’ ^^id U is just the 3x4 non-unitary mixing matrix 
in vacuum. The matrix U, that describes the effective neutrino mixing in matter, can be expressed 
as U = UX* which is also a non-unitary 3x4 matrix. 
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We find that it is easier to perform the diagonalization of TL' = UqUi/q , with 



(A5) 


by doing the unitary transformation T-L' = X'EX'^, and therefore the effective neutrino mixing 
matrix in matter is 


u = u{u^x'y = uulx'* 

^ “ 2 '^15 + ^le) 

~ -S14S24 + S15S25 + ^16'526 

^^“^ 14^34 + ^ 15®35 + ^ 16'536 


0 

1 “ 2 + '^25 + ■® 26 ) 

-S24'S34 + 'S25®35 + '526^36 ^ 


0 

0 

2 (^34 + 'S 35 + '^se) 



’24 


X'* 



(A6) 


In this new basis, we can write down the series expansion of the Hamiltonian as 

+7^'^^^ + ... , (A7) 

where and can be found in Eq. (A8) and Eq (A9), respectively. For the eigenvalues and 
eigenvectors, we also write as E\ = E^'^ + E^'^ + ... and + ... (for f = 1, 2 ,3,4), and 

the unitary matrix is defined as X' = {vi,V 2 ,v^,v^). 


22 



0 




= • 1 + 


Am|i 


+^NC 


5oq5i qCi q 6 


C23'®13'^13^ 


’23''23''13 


23 '’ 13 '' 13 ' 


^CC ^NC 


,2 2 

*12'-13 


Amii 


(^12^23 “ 'Si2'S23'5i3® *^) '*12^13 
“ (‘^12'523 + ^ 12 ^ 23 ^ 13 ^ **^) '5l2'^13 


0 Amli/Ar 


1^12^23 •Sl2'S23-®13® ) '512^13 

42^3 + «12S23Sl3 

— ^ sin 20^2 sin 2^23•5^3 cos S 

i sin 2023 (« 12 Sl 3 -c? 2 ) 

— I sin 20 j^ 2 Si 3 (cos 2^23 cos 6 — i sin (5) 
0 


(c12'523 + ^ 12 ^ 23 ^ 13 ^''^) '512^13 

i sin 2023 (si2«13-c?2) 

— I sin 20 j^ 2 Si 3 (cos 2023 cos 5 + i sin (5) 

«2 „2 I „2 „2 „2 

^12'^23 ^ '^12^23'^13 

+i sin 20^2 sin 2023^13 cos 5 


^14 + ■Sis + Sl 6 

'Si4'S24 + 'Si5'S25 + Sie-S^e 

Su'Sk + SlS'SsS + *16'*36 

— S 44 

14^24 + Sl5'S25 + Sl6'S26 

■*24 + '*25 + *26 

*24'*34 + ■*25'*35 + ■*26'*36 

“ ^*24 

I 4 S 34 + Sl5'S35 + Sl6'S36 

■*24*34 + ■*25*35 + ■*26*36 

*34 + *35 + *i 6 

“ ^*34 

£* 

*14 

£* 

*24 

£* 

*34 

- (Sl4 + *24 


+ycc 


- (si4 + ^15 + ^le) 0 0 ®14 

0 0 0 0 

0 0 0 0 

^14 0 0 5^4 



can be easily diagonalized. The resulting eigenvalues and the eigenvectors of are 


= E,- V^c , 

= E.-V^c + lVcc + l^i-^-C) 

-^ 3 °^ = -®1 “ ^NC + 'j^cc + 2 2E^ 


2E„ 


where C = y^l — 2AfjQ cos 20 ^ 3 , and 


(AlO) 


A'(o) = 


( 0 

Cl3 


S 13 C*' 

0 

1 

fi 

to 

CO 

“•^23^13® 

-iS 

'523'®13 

0 

•S 23 

~ 023^130 

-iS 

^23^13 

0 

V 0 

0 


0 

1 


where 


S ]^3 = sin 01 d = 


13 


Ci 3 = cos 013 = 


^ cos 2013 ^cc 

2 V C 


!}_ (^ I cos 2013 ^cc 
2 V C 


(All) 


Furthermore, the first order corrections to the eigenvalues and eigenvectors are given by 




= 


E 


w 


( 1 ) 

ji 


^ i ^ j 


(0) 


(A 12 ) 

(A13) 

(A14) 

(A15) 


where Inserting Eqs. (A9) - (All) into Eqs. (A14) and (A15), we can obtain 

the analytical expressions of E^^'^ 

^ = O' p,^^Ci2 + (^22C23 + ^33'523 “ -^23 sin2023) , 


Ek^> = 


1 Ami 


2 2E, 


'31 g2^ ( 2 + ^ ^CC COS 2013 


+ {yNC ~ ycc) ^llCl3 + VnC (^ 6^13 — sin 2013 ^ , 

Ef’ = (l - + (VVc - Vccmil, + live (li4 + , 

4 ^ = — — Vec) '^14 “ yNC ('^24 + 534 ) 1 (A16) 
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as well as the expression of 


= a sin 26* 


12 


0 


0 


(1 + A(j(j) S23'®13® * (^13'®13 + C13C13) C23 (Ci3'Si3 — ^13^3) 0336* 




^CC + (1 “ C') 


^CC 

(1 + Aq(j) 0235336 (^13533 + C33C33) S23 


^CC + (1 + C') 

^13^13 Sty-xC. 


i5 


13H31 •’23'= 


9 4 (-2 
0 

/ 


^CC + (1 “ C') 
0 


0 


sin 26*336 
2 C 2 


0 


’13 


-Ci36*'5 0 


\ 


0 •®23^13^ '^23^13 ii 

0 ^236336 C23S33 0 


^CC + (1 + C') 

0 oy 


(A17) 


/ 


\ 0 0 0 0 

Since the light sterile neutrino mass is supposed to be around the eV scale, we have safely omitted 
the terms proportional to sAm 33 /Am 43 . The effective neutrino mass-squared differences Am 23 , 


Am| 3 , Am |2 can then be written as: 


Amli = Am33 


2 / ^cc + (1 “ C*) 


— a 


„2 I „2 ^cc ^^13 (1 + 

C 12 + S 12 ^ 


+ i^NC ~ ^Cc) '^11^13 


-A 


NC 


^22 (^23 “ 'S23S13) + ^33 ('523 “ CI3S33) — R22, sin 26*23 (^ + 5 l 3 ) + S™ 2^33 | , 


Ami3 = + + 


— a 


C 12 s 


2 ^cc ^-'^s 2^33 (1 C) 


12 " 


2C 


+ i^NC ~ ^cc) 7^11513 


-A 


NC 


Amla = Am33 


722(523 CI 3 S 33 )-|-T 33 (c 23 S 23513 ) 7i23 sin 26*23 (1 + C 33 ) -Rfe sin 2033 } > 


9 1 COS 2^1 Q / . i \ r-r-i ^ 7\ 

C — 0:532 — {^NC ~ ^cc) 7^11 cos 2033 


+Ajq(j I Tfj cos 2033 + ^7?fe sin 2033 


(A18) 


and the effective mixing matrix U in matter can then be calculated by Eq. (A 6 ). Note that in this 
paper we do not order the eigenvalues of T-L according to their magnitude and the mass spectrum, 
since the ordering does not change the neutrino oscillation probabilities. 


Appendix B: Neutrino oscillation probabilities in matter 


In terms of effective neutrino masses and mixing in matter, the neutrino oscillation probabilities 
in matter j^) have the identical form as that in vacuum; 


6 /-) 
P{ ^ a- 




^ / 

fl* Ur, 

(e,= 

=1,2,3,4 (X]j=l,2,3,4 \ 

/ ^ ^ai^(3i 

i=l,2,3,4 


25 







\ 

4^Re 

'Ua^UpJKJUh 

sin^ Ajj ± 2 E Ini 

UarUpjKMi 

sin 2A^-, 1 

j>i 


j>i 


J 


(Bl) 


where Aj^ = Aruj^L/AE^ with Am|j = ?fi| — mf = 2 E^{E- — E-), and the lower sign is for the 
antineutrino oscillation probabilities. It should be noted that the mass indices i, j run over only the 
light neutrinos which can be kinematically produced in low energy neutrino oscillation experiments. 
For the Am|^-dominated accelerator neutrino oscillation experiments, oscillatory behaviors of the 
terms driven by Am|j (for i = 1,2,3) are averaged out. In this case the oscillation probabilities 
can be approximately written as 
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■21 


31 


32 


By inserting the results obtained in Appendix A, we are able to write down the neutrino oscillation 
probabilities —)• r'^) to the first order in both a = Am^i/Am\i and (for i = 1, 2, 3 

and j = 4, 5,6); 


p ~ 1 _ I _ 
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sin(C'A33) cos(C'A33) 


(B 4 ) 
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1 - 


1 2 (Sgg + Sgg) 


2s| 4 + sin 20 


23 


■ sin 2023 (1 “ 2T33) 


COS 2^ 

--P- (q;(1 + sin 2042 Si 3 cos 5 — (R^ sin 20^3 + 2 T^ (^CC + ^13)) ^Nc) 

^CC ^13 

+ 2 i ?23 cos 2023] [1 — COs(l + A(jq)A^i COs(C'A34)] 

C23 (1 “ 2T33) (2S2S cos 2043 sin(l + A34 + C23 sin^ 20^3 sin(C'A34)^ 


+ 


a sin 2042 2023S43 cos 6 

1 A.^^ cos 2043 


1 - 


cos 20, 


c 


(^CC + 251324^0 + ■S13) 


13 


sin(l + A(j(j)A2^i 


-2c; 


23 ' 


c 


sin(C'A34' 


A. 


— 2 Q;sf 2 C 23 ('^23 sin(l + A 34 — 20^3 cos 2043 sin(C'A 34 )^ 

—4i?43C23 sin 2043 ('®23 sin(l + ^C'c)'^31 “ C 23 cos 2043 sin(C'A 34 )^ 

+ 2 i ?23 sin 2023 (cos 2023 cos 2043 sin(l + Aqu)A^i — C 23 sin^ 2043 sin 


y4. fjr'^ Sin 2^90 sin 2^i q/tq /. 9\n./-. a \a 

+Ra - - - 2 - ([ci3 “ COS 2023 (^CC + ^13)] sin(l + A(^^)A. 


■31 


- 2 ci 3 C? 3 ^-^^ sin(C'A34)^ 


Ajyc sin 2023 

~^-^a ^ A 2 




"’"2 — COS 2023 ( 2 'Si 3 + ^CC ~ cos 2043 + ^^cc) 


^ sin^ 2043 


• sin(l + Aqq)A^^ — C23 sin^ 2043-^ sin(C'A34)^ 

- 4 c 23 sin 2043 


cos 2043 - 


- 2 ^bSin 2043 j ^ + -T 44 


^NC I ^NC ^CC 

C -J44 sin/c^43 — 


■ ^^23 sin(l + A(j(j)A^i — c^g cos 2043 sin(C'A34)^^ 

^Si2 ^^43 + - sin^ 29 i^A(j(^ — 042^ 
1 \ 


sin(CA34) 


+A31 sin^ 2023 


+ \ sin 2043 cos 2043 “ 


~ 2 ^b sin^ 2043 + ^NC 


+ ^^11 sin^ 2043 {^NC ~ ^cc) 


sin(l + A(j(j)A-^i cos(C'A34) 
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-Agisin^ 26»23 


(X I 5105 


2 „2 ^ + ^CC , ^2 


12'’13' 


c 


C^2 COS 2^]^^ 


“ 2 (COS ) AjYfj 


cos(l + Aq(j)/S.^i sin(C'A33) 


2^31 *^23 2^33 


as 


12 " 


1 cos 2^33 


c 


— ( i?f, sin 26*33 + ^^13 ) ^NC 


+ 2^11 ^°® 26*13 {^NC ~ ^cc) 


sin(C'A 3 i) cos(C'A 3 i) 


(B 5 ) 


P 




1 “ (^15 + ®16 + ®25 + ^ie) 


2 I 2 \ {l®15®25 + 'Sl6'526p + 


* " '^ ' ''^23 2013 (1 — Til — r 22 ) sin(C'A 3 i) 


+2as\2sl^ sin^ 20i3 ( cos 2^13sin(C'A3i) - A3i^—cos(CA3i) 


+ - 


a sin 2012 s™ 2023S13 


CA, 


cc 


cos 5 


/ -1 * \ A 1 cos 2010 . /^a a 

sin(l + a 4 ^c )^31 - ^ -sin(CA3i) 


-sin( 5 [cos(l + a4c-c)A3i - cos(C'A3i)]) 

p2s23 sin 2013 ^-^12 20^3 sin^(i 7 A31) -\- I12 cos(C'A3i)^ 


■A /\rp~* Sin 2^90 , 

' sin 20 


A, 


cc 


13 


Ra { sin(l + a 4 ^c )^31 - ^ sm(CA3i; 


+Ia (C0S(1 + A(jq)/S.^i — COs(C'A3i))] 
2 s 2 


+- 


’13 


c 


Ta (^sin(l + a 4 c.c-)A 3 i - ^ ^ sin(C'A 3 i; 


“-^23 (cOs(l + A(j(j)/S.^i — COs(C'A 3 i))]) 


+S 23 sin 20 i 3 ( — sin 40 


13 


Ji ^NC _ Ji ^NC ^CC 


c 


Ol" 


C 


sin(CA 3 i) - 2 '^A 3 i cos(CA 3 i) 




cos^ 2013 siii(C'A3i) + -C'A3 i sin^ 20i3 cos(C'A3i) 


sin(C'A3i) \ , (B6) 


P., 


1 


2 —3 rT2 TTTT 1 I'’15^35 t '’16'’36 


1 “ (^15 + ®16 + ®35 + ^ie) 


'STkS 33 + stfiSsfiP + 


C23 sin 2013 (1 “ ^11 “ ^22) siii(C'A3i) 


+2as?2ci3 sin^ 20i3 ( cos20i3^^ sin(C'A3i) - A3i^^— cos{C 


asin20in sin209oSio , 

Ia! AIS Id / ^ 


Ca 4 , 


cc 


/ t \ \ A 1 COS 2010 . /^A A 

sin(l + A(jc)^ 2 .i - ^ -sin(C'A 3 i) 


— sin (5 [cos(l + A(j(j)/S.y^^ — cos(C'A3i)]) 
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p, 


flT 


“1“2c 23 sin | -^13 cos sin({[/i^3]^) -\- /i3 cos(C'A3i)^ + i ?23 sin 2023 sin^ 20^3 sin(C'A33) 

I-A, 


+ 


./I ;\r^ Sin 2^00 . 

^^ ' sin 20 


A, 


cc 


13 


Ra ' sin(l + Aq(^)/S.^i — 


^cc 


C 


sin(C'A 3 i) 


+Ia (C 0 S (1 + A(j(j)/S..^i — COs(C'A 33 ))] 


+ 


‘ 2 ‘A 

^■^13 

c 


Ta ( sin(l + A^(j)l^^i sin(C'A3]^) 


~h3 (cos(l + A(j(j)/S.^i — cos(C'A 3;^))]) 
+C23 sin 20^3 sin 40;^3 
A 


jn Aj^(j ^ ^ Aj^(j A(J^ 


C 


11 " 


c 


sin(CA 3 i) - -CA 3 iCos(CA 3 i) 


+4i?i, 


^NC 

C 


1 


cos^ 2033 sin(C'A33) H—C'A33 sin^ 20^3 cos(C'A33) 


sin(CA 33 ; 


(B7) 


1 


1 “ ('®25 + '^26 + ®35 + ^ie) 

1 


I |'S 24 'S 34 C 23 “ (' 525^35 + 'S 26 ® 36 ) ^23 


2 |2 


+ 2 2023 (1 “ 7^22 — ^33) [1 — COs(l + ^(^(^)A33 COs(C'A33) 


+ COS 2033 sin(l + ^(^(2;)A33 sin(C'A33) — ^ 2033 sin (C'A33) 


--as ?2 sin 2 2023 sm^ 2033-g^ 


sin^l -j- 24 f^^)A 33 P cos 2033 sin(([ 7 A33) sin((i7A33) 


+ 


a sin 2032 2023 2023S33 cos S 


A, 


cc 


2 (1 — COs(l + A(j(j)/S.^i COs(C'A33)) 


"13 


+- 


^ CC + 2Si3^CC + ®13 

Ccfs 


sin(l + sin(C'A33) + 


1 — Ar<r' COS 20 


CC ^'^13 ^-2 


C 2 


sin (CA33) 


a sin 2019 sin 2090 5,0 sin (5 , , , , 

[cos(l + A(j(^)A^i — cos(( 7A33)] sin((7A33) 


CA 


cc 


— sin 2023 2033 


(^sin(l + A33 + cos 2033 sin(C'A33)^ 

—la (cOs(l + ^(7 c')A33 — COs(C'A 33))] sin(C'A33) 


+i?23 sin 2023 [cos 2023 (1 ~ COs(l + A(jq)A'^i COs(C'A 33 )) 


+ cos 2023 cos 2033 sin(l + Aq(j)A-^i sin(C'A33) + 5^3 sin^ 2033 sin^(C'A33) 


“t“ Sin ‘ 2 i 0 o‘il‘ 


23-'23 


sin(l cos((^z^3]^) “h cos cos(l ^(27C')^3i sin(C^z!\3^ 
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,2 „2 ^CC + 2 Si 3 ^CC + 1 „2 


'Sl 2 Sl 3 ‘ 


C 2 


— Ci2 sin(l + v1(^C')A3i cos(C'A3;^) 


. 2 2 
“ I •Si 2 'Si 3 


1 + A 


CC , ^2 


c 


c^2 COS 2(9i3 cos( 1 + ylc-t;)A 33 sin(C'A 3 i) 


+ 2^12 2 ( 9 i 3 


1 COS 2^1 o / K \ / ^ K \ 

- ^ -sin(CA 3 i) cos(CA 3 i) 


+ l^^ -4^csm4A3Sin2g,3 ^ A^^)A33 cos(CA3,)) 


1 


“ {^CC + ^13) sin(l + Auq)/S.^i sin(C'A33) — ®™^(C'A3i) 


+A,, [{^cc + ■513) (1 “ cos(l + Aqq)A^i cos{CA^i)) 


A A 


2^33 ~t“ Aq^ — 3A^q cos 2^^3 -j- ‘2Aqq 

2 C 


sin(l + Aq<q)A^i sin(C'A3i) 


-- sin^ 2013 (1 + ^cc) sin^(<^^3i) 


yljvcsin2023^sin2013 ) [cos(l + Acc)^3i “ cos(CA3i)] sin(CA3i) 


-Rb COS 2013 - 2 ^f'®™ 2 '^i 3 


^ + sin 2^,3 -lee ) 3 i „2 2^^^ 2033 


sin(l + A(^q)A^i + cos 20i3 sin(C'A3i) sin(C'A3i) 


+ o^3i^fe2liv^sin2 2023 sin2013 cos 20i3 sm(l + ylc'^)A3i cos(C'A3i) 


+ cos(l + A(j(^)A^i sin(C'A3i) — sin^ 20i3 sin(C'A3i) cos(C'A3i) 


“^^31 [^b^NC ~ ^11 i^NC ~ ^cc)] sin 2023 s™ 20i3 


sin( 1 -\- “h cos 2^]^2 sin(^{[/ “^31) cos) 


+ 2 ^ 3 iA; 24 - 7 VC sin^ 2023 [sin(l + A(j(^)A^^ cos(C'A3i) 
“ 1 “ cos 2013 cos(l “ 1 “ AfjQ^jA^-^ sin^(^A3i) ^ . 


(B 8 ) 


The above analytical approximations become invalid at the relatively large values of LjE^ or for 
the Am|i-dominated oscillation. Note that, one cannot obtain the vacuum oscillation probabilities 
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Pi^u 1^) from in Eqs. (B 3 ) - (BIO) by directly setting Vq(j = = 0 , because 

the expansion of % in Eqs. (A 4 ) - (A6) would be not appropriate in the limit of vanishing V(j(j 
and VjqQ. In addition, one can obtain the corresponding antineutrino oscillation probabilities 
= Pii'a using the relation 

~^CC’ ~^Nc) = Papi^-: ^Nc) ■ (^ 9 ) 

Because of the existence of unitarity violation, the “to all” probabilities Pii'a —^ ^e,^i,T) = ^ae + 
P^^ + Par a = e, n,T are in generally not the unity, but displayed as 

P{^e ~ 1 “ 2 s ?4 + Sm 20^3 (5^4 - 5145^3 - 5140^3) sin 2043 + 2 -Rfe COS 2^43 

+i?23 sin 2023 sin2043 sin^(C'A34) +/^ sin2043 sin(2C'A34) , (BIO) 



1 - 25^4 - (s?4 - - 5^40^3) 5^3 sin^ 2043 sm2(C'A34) 

sin^ 2023 { 4 a - 4 a) [1 - cos(l + Acc)^zi 


“h COS ^ 20 -^^ sin^l -|- 


- sin 2043 Pa sin 2023 sin(l + A^c-) A34 + 2 i?f,S 23 cos 2043 sin(CA34) 
—/„ sin 2023 cos(l + ^cc)A 3 i + 24^23 cos(C'A 34)] sin(C'A34) 


+i ?23 sin 2023 [cos 2023 (1 ~ COs(l + COs(C'A34)) 

~t“ cos 2023 cos 2043 sin(l -1- A^^)A34 sin(( 17 A 34 ) -t- S23 sin 2043 sin ((I7A34) 

+/23 sin 2023 [sin(l + A,^,^)A34 cos(C'A34) 

+ cos 2043 cos(l + A(^(^)A34 sin(C'A34) ) (Bll) 


Pifr ~ 1 — 2 s |4 — (5^4 — S24S23 “ SI4C23) C23 sin^ 2043 sin^(C'A 34 ) 

— - sin^ 2023 (^24 “ •SI 4 ) [1 — COs(l + A(jq)A^i COs(C'A 34 ) 

-\~ cos 20 ^^ sin^l -|- sin^C^zX^i) 

+ sin 2043 Pa sin 2023 sin(l + A^c-) A34 - 2 R^cl^ cos 2043 sin(C'A34) 
—4 sin 2023 cos(l + ^cc)A 3 i — 24C23 cos(C'A 34 )] sin(C'A 34 ) 
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— i ?23 2^23 2^23 ~ COs(l + Aq(j)A^i COs(C'A3;^)) 

+ COS 2^23 COS 2 ( 9^3 sin(l + AfjQ)A.^^ sin(CA 3 ;^) — (2 + CI 3 ) sin^ 2^13 sin 2 (C'A 3 i) 


—/23 sin 2023 [sin(l + A 33 cos(C'A 33 ) 


“h COS “ 20 ^^ COs(l “h -^(7c)^31 


(B12) 


where we have used the relation 


4 ) = Pcai^^ 4 ) > 




0/3 ^ij) 


(B13) 


in matter of constant density. One can find that these summed probabilities —>■ T^e,iJ,,T) 

independent of the neutral-current potential Vj^^. 

Considering the smallness of S 33 , one can make further simplifications for these oscillation 
probabilities. To the order of S 33 , we have 


C 


|1 “ ■^cc\ + 


2A. 


\l- A 


■813 — 7 


cc\ 


(1 “ ^cc) + 


„2 

^ ^cc 


where 


7 = sign(l-^cc') = ^ 


1-A 


cc\ 


Then we have 

sin(CA3i) PS 
COs(C'A3i) Ri COs(l — A(^(^)A33 — 2A33S33 


sin(l — A(j(j)A.^^ -|- 2 A 33 S ^3 


^cc 


1 - A, 


■ cos(l — Aq(j)A^^ 


cc 


A 


cc 


1 - A 


■ sin(l — Aqu)A^i . 


cc 


(B14) 


(B15) 


(B16) 


If we neglect the double-suppressed terms 0{s\2^a) or 4 ) or 0 ( 5 ^ 3 ), the more concise ap¬ 

proximate formulas for the neutrino oscillation probabilities in Eqs. (11) - (16) and Eqs. (17) - (19) 
can be obtained. 
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